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ABSTRACT 

The local axisymmetric stability of hydro dynamical and magnetized, nearly- 
Keplerian gaseous accretion disks around non-rotating black holes is examined in the 
vicinity of the classical marginally-stable orbit (at radii ^ Rms)- An approximate 
Paczynski-Wiita pseudo-Newtonian potential is used. Hydrodynamical disks are 
linearly unstable inside a radius which differs slightly from the classical -Rms value 
because of finite pressure and radial stratification effects. Linear stresses associated 
with unstable hydrodynamical modes vanish exactly at the radius of marginal stability 
and are generally positive inside of that radius. When a magnetic field is introduced, 
however, the concept of radius of marginal stability becomes largely irrelevant because 
there are linearly unstable magneto-rotational modes everywhere. Associated linear 
stresses are positive and continuous across the region of hydrodynamical marginal 
stability, even for large-scale "hydro- like" modes subject only to weak magnetic tension. 
This conclusion is valid for arbitrarily thin disks (in ideal MHD) and it does not 
require a large-scale "radially-connecting" magnetic field. Results on hydrodynamical 
diskoseismic modes trapped in deep relativistic potential wells should be revised to 
account for the short-wavelength Alfven-like behavior of inertio-gravity waves in 
magnetized disks. 

Subject headings: accretion disks — hydrodynamics — MHD — instabilities — 
turbulence — black holes 

1. Introduction 

Over the past few years, the question of the nature of accretion in the vicinity of black 
holes has received considerably renewed attention. This interest stems from a combination of 
new observational data and a deeper theoretical understanding of the accretion problem. On the 
observational side, significant progress has been made since a new generation of powerful X-ray 
telescopes (Chandra and XMM-Newton) became available. These instruments have been used to 
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strengthen earlier claims of evidence for general relativistic effects arising deep in the potential 
well of spinning black holes, in the form of broad iron emission line profiles (Iwasawa et al. 1996; 
Dabrowsky et al. 1997; Branduardi-Raymont et al. 2001; see also Reynolds &: Begelman 1997; 
Young, Ross & Fabian 1998). In addition, it has been suggested that in some systems one actually 
witnesses black hole spin energy being tapped to power part of the observed X-ray emission 
(Wilms et al. 2001, Miller et al. 2002). 

On the theoretical side, efforts to reconsider classical results (see, e.g., Novikov k, Thorne 
1973; Page & Thorne 1974, Shapiro & Teukolsky 1983) on the marginally-stable orbit and the 
zero-torque boundary condition at i^ms have been motivated by the compelling case made for 
magnetic fields and MHD turbulence being the driving mechanism behind accretion (Balbus &; 
Haley 1991; 1998; Balbus 2003; Blaes 2003). Recent theoretical efforts on near-black-hole accretion 
have focused on analytical models of steady-state accretion (Krolik 1999; Gammie 1999; Agol 
k, Krolik 2000; Paczynski 2000; Afshordi & Paczynski 2003) or on global numerical simulations 
of magnetized flows in an approximate pseudo-Newtonian potential for non-rotating black holes 
(Hawley & Krohk 2001; 2002; Reynolds & Armitage 2001, Armitage, Reynolds & Chiang 2001; 
Krolik & Hawley 2002). 

Gammie & Popham (1998; see their footnote 7) and Araya-Gochez (2002) have studied some 
important aspects of the linear stability problem for magnetized disks around black holes in full 
general relativity, but it seems that there has been no attempt to use this approach to address 
some of the debated points in the literature on the subject of near-black-hole accretion (e.g. 
validity of zero-torque condition, role of disk thickness, differences between hydrodynamical and 
magnetized disks) . This is undoubtedly because there are issues that a local, linear analysis cannot 
fully address and which require global numerical simulations. Nonetheless, as I argue below, such 
an analysis still provides valuable insights on some of the key aspects of the accretion problem in 
the vicinity of a black hole. In §2, I describe the formalism used for the linear analysis. I present a 
comparison of results for hydrodynamical and magnetized thin accretion disks in §3. Implications 
and limitations of this work are discussed in §4. 



2. Formalism 

We work in cylindrical coordinates {R,(j),Z). Without loss of generality, we restrict our 
analysis to a basic state disk with a purely vertical magnetic field, B = (0, 0, B), and a rotation law, 
ft = (0, r2(i?), 0), that is constant on cylinders. We consider axisymmetric Eulerian perturbations 
(denoted by a prefix 6) with WKB space-time dependence e*('^^~'^*). We work in ideal MHD (i.e. 
microscopic diffusion processes are ignored) and we neglect weak circulations of the basic state, 
such as that caused by a finite fiuid viscosity. 

The formalism used for this analysis is essentially the same as the one developed by Balbus 
&; Hawley (1991; 2002). Standard MHD equations and expressions for the leading-order WKB 
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terms (within the Boussinesq approximation) can be found in the first of these two references. 
Two of these equations for hnear perturbations are useful for the discussion below and they are 
reproduced here: 



-ludvs + dVR — - ikB—^ 
2il Airp 



0, 



(1) 



-iuSBji — ikBSvR = 0. 



(2) 



Linear perturbations satisfy the dispersion relation 



(3) 



where 



VA = B/^/4^^p, u^=u?-{kvAf, 



(4) 



and p is the mass density of the fluid. Unstable modes are best characterized by their 
growth rate, 7 = —iuj. The Brunt- Vaisala frequency and the epicyclic frequency appearing in the 
dispersion relation are defined, respectively, by 



N^ = -- ^ 



«2 = 
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where P is the gas pressure and 5/3 is the adiabatic index of the gas. 
Following Balbus & Hawley (2002), the growth rate, 7, of unstable modes obeys 
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The wavenumber (k) and growth rate of the fastest growing mode are given by 
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The R(f) component of the linear perturbation stress tensor, which measures the angular 
momentum flux associated with unstable modes (Balbus & Hawley 2002; 1998), is 
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TR(j, = p [SvrSv^ - SvarSvacpI , 
where the first term is the Reynolds stress and the second term is the Maxwell stress: 
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and D = 1 + (kvA/'j)'^- In what follows, we will focus on the dimensionless stress 
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While Balbus & Hawlcy (2002; see also Narayan et al. 2002) were interested in the general 
stability and transport properties of accretion flows with strongly destabilizing radial stratifications 
{N^ ~ —O!^), our specific interest here is in the possible change of stability properties caused by the 
gravitational potential of a non-rotating black hole in the vicinity of the classical marginally stable 
orbit (-Rms)- We approximate the Schwarzschild metric with the pseudo- Newtonian potential of 
Paczynski &: Wiita (1980), which implies a Keplerian rotation law 



GMbk 
R{R - RgY ' 



(12) 



where G is the gravitational constant, Mbh is the central black hole mass and Ry = 2GMbh/c^ 
is the Schwarzschild radius of the non-rotating hole (c is the speed of light). This approximate 
potential is known to reproduce accurately the location of the classical marginally-stable 
{Rms = 3-Rg) and marginally-bound (-Rmb = 2i?g) orbits. We focus our attention on the stability 
of nearly-Keplerian gaseous accretion disks, i.e. gaseous disks in near rotational equilibrium. In 
the absence of any accretion (i.e. infall), the equilibrium (basic state) rotation law of such a disk 
is governed by the radial component of the Euler equation. 



RQ^ 



pdR 



(13) 



The assumption of near Keplerianity means that radial pressure gradients are small compared 
to the two other terms in the above equation: O ~ . This is equivalent to requiring that the 
disk geometrical thickness be small compared to its radius {H/R ^ 1). A further consequence 
of this assumption is that, absent any variations of the density or temperature over a Icngthscale 
comparable to the disk pressure scale height, H, the epicyclic frequency, k^, is close to the 
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Keplerian value as well. In what follows, we will neglect any such deviation of from local 
Keplerianity. 

Since our basic state disk configuration is, by assumption, non-accreting, it would seem 
natural to assume that any radial stratification (that would in general arise from accretion 
and associated dissipation) is absent. It is possible, however, to imagine situations in which a 
non-accreting gaseous disk still possesses a non-zero radial stratification, as could be the case 
for instance in the presence of external sources of irradiation. For this reason, we generalize 
our discussion to the cases with: (i) no radial stratification (AT^ = 0), (ii) a stabilizing radial 
stratification {N"^ = O.lJl^) and (iii) a destabilizing radial stratification {N^ = —O.lfP). The 
values ibO.iri^ are intentionally chosen to be large (for a thin disk), to better illustrate the effects 
of finite N"^ values on stability. Reducing these values to more realistic levels would not affect the 
results qualitatively, but it would make the identification of the effects on figures more difficult. 

3. Results 

3.1. Stability of Hydrodynamical Disks 

While the formalism described in the previous section has been developed for magnetized 
disks, the hydrodynamical limit can be recovered by enforcing kvA = 0. The dispersion 
relation reduces to second order and the growth rate of unstable modes becomes independent of 
wavenumber (k): 7^ = —(N'^ + ac^). One recognizes one of the classical Solberg-H0iland criteria 
for axisymmetric stability (e.g. Tassoul 1978) emerging from this expression for 7^. Around a 
black hole, the Keplerian rotation law becomes unstable inside of the classical marginally-stable 
orbit, i?ms = 3i?g, where < (for all wavcnumbcrs k). While this is an exact result for a test 
particle (or a dust disk), hydrodynamical disks behave somewhat differently. First, the value 
of the disk epicyclic frequency, k^, is slightly modified by the presence a finite radial pressure 
gradient (eq. [13]), a small effect that we have chosen to ignore for simplicity. Second, the presence 
of radial stratification modifies the radius of marginal stability (i.e. the radius at which 7^ = 0). 
Marginal stability occurs at a radius somewhat larger (smaller) than the standard Rms = ^Rg in 
the presence of a destabilizing (stabilizing) radial stratification. These modifications remain small 
for cold/thin disks (with N"^ < n'^). 

Dimensionless Reynolds stresses associated with linearly unstable hydrodynamical modes 
take the (wavenumber independent) form: T^^ = —{k?/2Q'^) x (^2/7). In the absence of radial 
stratification, 7 = \/— k^, the dimensionless stress is everywhere positive within the radius of 
marginal stability and it vanishes exactly at that radius. When a finite radial stratification is 
present, however, the degeneracy between the growth rate, 7, and the epicyclic frequency, k^, is 
broken, and T^^ diverges at marginal stability (7 = but 7^ 0). The divergence is positive 
for a stabilizing radial stratification and negative for a destabilizing one. Therefore, in a disk 
with A''^ < 0, there is a small region with negative stresses right inside of the radius of marginal 
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stability. These results are graphically illustrated in figure 3 (see below). 

The divergence of the dimensionless stress does not apply to the physically-relevant dimensional 
stress, however. Equation (1) for linear perturbations shows that, in a hydrodynamical disk (with 
B = 0), when 7 = — icj ^ while 7^ 0, 5vji must ^ as well. [This requirement disappears 
if = at marginal stability (i.e. absent any radial stratification), but there is no divergence 
of in that case because 5vf{ can remain finite (see also Fig. 3).] The fact that 6vr ^7^0 
at marginal stability is precisely what maJces the dimensionless stress, T^^, diverge (if iV^ ^ 0). 
On the other hand, the dimensional stress, Tji^ oc x {Svr)^, vanishes exactly at marginal 
stability, where 6vr — 7 — 0. All this guarantees that the physically relevant linear stress always 
vanishes at marginal stability, whether there is radial stratification or not. Note that a similar 
situation occurs at marginal stability {Sv^i ^ 7 ^ 0) in a magnetized disk, as equation (2) for 
linear perturbations shows. Despite the possibility of divergence at marginal stability, we found 
that using T^i remained very useful because of its dimensionless nature.^ 



3.2. Stability of Magnetized Disks 

Figure 1 shows growth rates of unstable modes, as a function of dimensionless wavenumber, 
kvA/^, at four different locations in the vicinity of the classical marginally stable orbit, i^ms = 3Rg. 
Only the case without radial stratification (iV^ = 0) is shown because it was found that a value 
of A'^^ << only slightly reduces or increases growth rates in an approximately wavenumber 
independent way. For each of the curves in figure 1, the mode with maximal growth is easily 
identified at kvA ~ The growth rate of this most unstable mode increases as one approaches the 
black hole horizon because the shear and the epicyclic frequency become ever more destabilizing. 
Values of the growth rate for the most unstable mode shown in figure 1 are found to be in very 
good agreement with corresponding fully general relativistic calculations for Schwarzschild black 
holes (Araya-Gochez 2002). 

One observes a qualitative change in the behavior of small wavenumber (large scale) modes as 
a function of radius. At radii R > 3Rg (for N'^ = and a strictly Keplerian epicyclic frequency), 
the growth rate of these modes tends to zero in the limit kv^ 0. Within 3Rg, however, the 
growth rate tends to a finite value, which is 7 = \/ —n^ . This is exactly the growth rate of 
unstable modes in a hydrodynamical disk at the same location within the radius of marginal 
stability. Therefore, large scale modes characterized by weak magnetic tension (small kvA) 
behave asymptotically like their hydrodynamical counterparts. Note that this is the only physical 
significance attached to the hydrodynamical radius of marginal stability in a magnetized disk. 
Because there always are unstable modes, on scales small enough for magnetic tension to be 



^Narayan et al. (2002) use the dimensionless quantity tr<t, = T^^ x (7/f2), which remains finite at marginal stability 
where the dimensional stress, Tr^, vanishes (and T^^ diverges). 
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important, there is no effective radius of marginal stability in a magnetized disk. Finite radial 

stratification or pressure effects can modify somewhat the exact radius at which small wavenumber 
modes become unstable (from the canonical 3Rg value; sec discussion of Fig. 3 below), but they 
do not qualitatively alter the stability or stress properties. 

Figure 2 shows the total (Reynolds + Maxwell) dimcnsionlcss stresses, TJj^, associated with 
linearly unstable modes, as a function of the dimcnsionlcss wavenumber, Uva/^, at the same 
four locations in the vicinity of the classical marginally stable orbit (same notation and N'^ = 
assumption as in Fig. 1). To illustrate the vanishing of the dimensional stress at marginal stability, 
we have multiplied the dimensionless stress by (in this figure only), since we know that 

5vR ^ 7 ^ at marginal stability (see §3.1). Again, a qualitative change in the stress across 
the hydro dynamical radius of marginal stability occurs only for large scale modes. The stress is 
asymptotically zero outside of 3i?c, and asymptotically finite inside of 3i?g, in the limit kvA 0. 

Figure 3 summarizes the picture that emerges for linear perturbation stresses in magnetized 
and hydrodynamical disks in the region around the classical marginally stable orbit. The stress 
for the fastest growing mode in a magnetized disk is positive and continuous throughout the 
region of interest. The introduction of a stabilizing or destabilizing radial stratification does 
not qualitatively change this result. On the contrary, the (wavenumber independent) linear 
perturbation stress in a hydrodynamical disk is positive only within the radius of marginal 
stability, whose location is affected by the introduction of radial stratification (see lower short- and 
long-dashed lines in Fig. 3). At marginal stability, hydrodynamical stresses vanish (the apparent 
singularities of T^^ in Fig. 3 do not apply to the dimensional stress, T^^^, as discussed in §3.1). 
The dotted line in figure 3 shows the linear perturbation stress for a representative "hydro-like" 
mode in a magnetized disk, i.e. a mode with [kvA/^)'^ = 10~^ which is subject only to weak 
magnetic tension. Even for this "hydro- like" mode, whose properties become very similar to those 
of purely hydrodynamical modes as i? ^ 2i?g, the linear perturbation stress remains strictly 
positive and continuous throughout the region of interest (as does the growth rate; see Fig. 1). 



3.3. Waves in Magnetized Disks 

While it is not the main focus of this paper, it is worth emphasizing here what can be learned 
from the stable branch of the dispersion relation (eq. [3]). Indeed, following early investigations 
by Kato &; Fukue (1980) and Okasaki, Kato &; Fukue (1987), it has been suggested that global 
hydrodynamical modes in unmagnetized accretion disks could be trapped in the vicinity of black 
holes (Novak & Wagoner 1991; 1992; 1993; Perez et al. 1997; Kato 1990; 1993; 2002; 2003; Li, 
Goodman & Narayan 2003). Mode trapping in this context relies on a purely relativistic effect: 
the epicyclic frequency reaches a maximum, before vanishing at the classical marginally-stable 
orbit, as one goes deeper and deeper into the black hole potential well (or that of a very compact 
neutron star for that matter; see, e.g.. Wagoner 1999 for a review). One of main motivations for 
elucidating the properties of such diskoseismic modes is that they may be at the origin of a variety 
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of quasi-periodic oscillations (QPOs) observed in compact accreting systems (see, e.g., van der 
Klis 2000 for a review on QPO observations). 

Given the restricted geometry considered, as well as the local and axisymmetric nature of 

the analysis, our dispersion relation (eq. [3]) describes only a very small subset of all diskoseismic 
modes discussed in the literature and global trapped modes cannot be studied explicitly. 
Nonetheless, the dispersion relation is interesting because it does describe, in the magnetized 
context, inertio- gravity waves^ which are often regarded as the most observationally relevant 
diskoseismic mode (Wagoner 1999). 

Figure 4a shows the frequency, w, of wave solutions to our dispersion relation, as a function 
of the dimensionless wavenumber, Uva/^- Hydrodynamical analogues to these waves can be 
identified by taking the limit kvA — ^ 0. Outside of the hydrodynamictd radius of marginal stability, 
we find that the asymptotic value of the wave frequency is w — > \//?~+~/V^ when kvA 0. Clearly, 
these waves correspond to the (incompressible) inertio-gravity waves of hydrodynamical theory.'^ 
As expected, the effect of introducing a finite radial stratification (not shown here) is to increase 
or decrease the wave frequency relative to the epicyclic frequency, (in addition to changing the 
location of marginal stability). However, it is obvious from figure 4 that the behavior of these 
waves is very diflFerent in the large-fcv^ limit, i.e. when magnetic tension becomes important. Since 
we are working in the Boussinesq approximation, fast magnetosonic waves are effectively filtered 
out from the dispersion relation. We also know that unstable modes in differentially-rotating disks 
are related to the slow MHD mode, so the waves described by figure 4 must correspond to the last 
of the three classes of MHD waves, namely Alfven waves (as expected from the larger frequencies 
of these waves relative to the frequencies of slow MHD waves - described by the other branch of 
our dispersion relation - in the limit Uva/^ 3> 1; see Balbus & Hawley 1998 for a discussion of 
MHD waves in the accretion disk context). 

Figure 4b illustrates two interesting properties of inertio-gravity- Alfven waves^ in the vicinity 
of the hydrodynamical radius of marginal stability. While inertio-gravity waves become purely 
evanescent within Rj^s in unmagnetized disks, the corresponding waves in a magnetized disk 
continue to exist within .Rms) albeit with very low frequencies (they become evanescent only 
asymptotically, i.e. in the limit kvA 0). Even if these waves do not become strictly evanescent 



^We adopt the standard nomenclature of geophysical fluid dynamics by referring to (internal) gravity waves in a 
rotating medium as inertio-gravity waves (see, e.g., Holton 1992). 

*Notc that because of the restricted geometry adopted (with only kz finite), the analogues of hydrodynamical 
diskoseismic waves with lu < k (absent any radial stratification) are not described by our dispersion relation. 

^Note the following interesting property of inertio-gravity modes in a magnetized medium. When they are stable, 
these modes are related to Alfven waves, as we have just discussed. However, when they are unstable, these modes 
switch branch in our dispersion relation and become then related to slow MHD modes. This can bo seen by comparing 
figures 1 and 4 or, equivalently, by comparing the two (stable and unstable) branches of the dispersion relation in 
the limit kvA — > 0, both inside and outside of the hydrodynamical radius of marginal stability. 
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(w — 0) inside of Rms, the long- wavelength behavior they exhibit is broadly consistent with the 

hydrodynamical limit in the sense that there is a maximum and a decrease of the wave frequency 
as one approaches -Rms- Note that, contrary to the purely hydrodynamical case, the location of 
the frequency maximum is not fixed at i? = ARg in magnetized disks but it becomes a function of 
wavenumber. Waves with different wavenumbers would then be trapped in different regions of the 
disk. 

All this suggests that some of the results on trapped inertio-gravity waves in hydrodynamical 
disks (due to the presence of a maximum of the epicyclic frequency, k^) may remain valid in 
magnetized disks, provided one is allowed to consider long enough wavelengths for kvA/^ to be 
<^ 1 {a condition which depends on the magnetic field strength explicitly). On the other hand, 
when magnetic tension becomes important (i.e. if Uva/^ is not small), inertio-gravity- Alfven 
waves behave much more like Alfven waves. Figure 4b shows that, in that regime, there is no 
significant change in the wave properties as one crosses the hydrodynamical radius of marginal 
stability. The frequency maximum disappears and, as a result, mode trapping is no longer 
expected. 

Our short digression on inertio-gravity- Alfven waves indicates that a more general treatment 
of waves in magnetized disks is much needed to address the issue of global mode trapping (including 
compressibility effects allowing the existence of fast magnetosonic waves). Our simple analysis 
suggests that results established in the hydrodynamical limit do not carry over automatically to 
magnetized disks, especially when magnetic tension becomes important. Another difficulty that 
becomes apparent when studying magnetized disks is that all the waves coexist with the unstable 
slow MHD (magneto-rotational) mode, so that ignoring the resulting turbulence in the treatment 
of the background disk may be a serious oversimplification. 

4. Discussion and Conclusion 

The classical description of thin accretion disks around black holes involves the action of an 
"anomalous viscosity" far from the hole, which (slowly) brings gas down to the marginally stable 
orbit, i?ms- At that point, particle orbits become unstable and the gas is expected to flow into 
the hole on essentially free-fall orbits. This rapid infall, whose origin is the dynamical instability 
at Rms, niust then cause the (steady-state) gas density within i?nis to be very small, so that 
fluid stresses should essentially vanish. Another reason why one expects stresses to disappear 
somewhere outside of the black hole event horizon is that the infall speed must reach the speed 
of light at the horizon according to general relativity and therefore the flow must lose dynamical 
contact with external regions once it crosses the sonic (or fast-magnetosonic) surface. 

Several elements of the classical picture should be revised now that we have a deeper 
understanding of the process responsible for accretion in the disk (i.e. the "anomalous viscosity" 
of the classical theory). Indeed, the magneto-rotational instability is a dynamical instability itself. 
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It relies on a rotational imbalance that is also what fundamentally leads to free-fall within the 

marginally stable orbit according to classical theory. Clearly, the fact that magnetized Keplerian 
accretion disks are subject to the same dynamical, magneto-rotational instability far away from 
the black hole and within the classical marginally stable orbit, and the realization that this 
instability is responsible for accretion far from the hole in the first place, suggest that the classical 
marginal stability theory may be inadequate to describe magnetized accretion disks. 

The linear calculations presented in this paper are an attempt to formulate these objections in 
a more rigorous way. There is no doubt that the classical stability theory is exact for test particles 
(or dust disks). Hydrodynamical calculations show, however, that positive linear perturbation 
stresses resulting from rotational imbalance (Rayleigh instability) in a nearly-Keplerian gaseous 
accretion disk are possible within (and only within) the hydrodynamical radius of marginal 
stability. This result is at odds with the classical description which postulates that the "anomalous 
viscosity" operates in the disk outside of i?ms but not inside of Rms- Even the hydrodynamical 
analysis may have little to do with astrophysical systems, however, since those are expected to 
contain magnetized disks. Once a magnetic field is introduced, black hole accretion disks become 
dynamically unstable to the same magneto-rotational modes everywhere. Disks are more unstable 
within the hydrodynamical radius of marginal stability because of the increasing shear and ever 
more destabilizing epicyclic frequency, but there appears to be no reason (at least in linear theory) 
why one should invoke a sharp transition at or about the classical marginally stable orbit, as 
illustrated in figure 3. 

It is interesting that this conclusion is independent of the disk geometrical thickness. Indeed, 
the only requirement for our analysis to be valid (within ideal MHD) is that the magnetic field be 
weak enough for the wavelength of at least one unstable mode to fit within the disk scaleheigth. 
This is roughly equivalent to the requirement of a subthermal magnetic field {va < cs; see, e.g., 
Balbus & Hawley 1998). In our analysis, one is therefore free to chose the disk sound speed (or 
equivalently disk thickness) as small as desired, as long as it remains finite. The results remain 
unchanged, provided the basic state magnetic field in the disk is made weak enough to allow for 
the presence of unstable modes. 

There have been arguments presented in the literature for and against the classical picture, 
and in particular related to the issue of vanishing stresses within i^ms (zero-torque boundary 
condition). Krolik (1999; see also Agol &; Krolik 2000) and Gammie (1999) have argued that 

stresses could remain large across R^g in a magnetized disk, while Paczynski (2000) and Afshordi 
&: Paczynski (2003) have argued in favor of the (near-) vanishing of stresses at radii R < Rms, 
provided the disk is geometrically thin (i.e. H/R <C 1). Qualitatively, our linear analysis indicates 
that stresses should not vanish exactly anywhere around i?ms in a thin, magnetized disk, and in 
that sense it supports the claims of Krolik (1999) and Gammie (1999). As a corollary, it suggests 
that the classical R^s location loses much of its significance in magnetized disks. However, to the 
extent that the zero-torque issue is a quantitative one (i.e. large vs. small stresses inside of i?ms)) 
it is unclear whether linear stability theory can provide a useful answer. One should also be aware 
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of important differences between our approach and the methods of previous authors. While these 
authors focused on steady-state models (and some of their conclusions rely on that property), our 
stability analysis is intrinsically time-dependent (a characteristic which may be viewed as a plus). 
On the other hand, we do not solve the global near-black-hole accretion problem, and in doing so, 
we neglect the important role of large infall velocities near the black hole event horizon {vr c) 
which determine the location where dynamical contact is lost (at the "sonic" surface) and where 
stresses should vanish. 

Another interesting property of our linear analysis is that it is local. The role of a large-scale 
magnetic field has at times been invoked to couple regions of the flow within R^s to regions more 
distant from the black hole and thus justify the idea that stresses may not disappear inside of Rms 
in magnetized disks (see, e.g., Krolik 1999; Agol & Krolik 2000). Our analysis, because it is local 
in nature, shows that such a large-scale "radially-connecting" magnetic field is not required for 
stresses to remain finite across the hydrodynamical radius of marginal stability. 

Clearly, there are a number of elements of the near-black-hole accretion problem that the 
linear analysis cannot capture. Most importantly, it is the global MHD turbulence resulting from 
the non-linear development of the magneto-rotational instability that is ultimately relevant for 
the observational properties of accreting black holes (see Krolik &; Hawley 2002 for a discussion 
of the issues in that regime). It is also significant that, while magnetized turbulent accretion 
disks have been shown to remain relatively close to rotational equilibrium (i.e. Kcplcrian) at large 
distances from the black hole (see, e.g., Hawley, Balbus &: Stone 2001), numerical simulations 
show (and the infall constraint vr ^ c guarantees) that accretion must be far from rotational 
equilibrium in the vicinity of the black hole event horizon (see, e.g., Hawley Sz Krolik 2001; 2002). 
This is an indication that linear perturbation results based on a nearly-Keplerian basic state must 
be interpreted with caution in that region. One should also keep in mind that our analysis was 
carried out for strictly axisymmetric perturbations and that non-axisymmetric perturbations may, 
in general, have different stability and transport properties. Nevertheless, our results suggest 
that there is no reason for stresses to vanish in the vicinity of the classical marginally-stable 
orbit, irrespective of the geometrical thickness of thin accretion disks and whether large-scale 
"radially-connecting" magnetic fields are present or not. 
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Fig. 1. — Dimensionless growth rate, 7/^2, of linearly unstable modes as a function of dimensionless 
wavenumber, kvA/Q, in a magnetized disk without radial stratification. The various curves 
correspond to different locations in the vicinity of the classical marginally-stable orbit at i?ms = 
(solid: 4:Rg; long-dashed: 3Rg; short-dashed: 2.5Rg; dotted: 2Rg). 
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Fig. 2. — Dimensionless stress, T^^ [x^'^/Q?), of linearly unstable modes as a function of 
dimensionless wavenumber, kvA/^, in a magnetized disk without radial stratification. Notation 
is identical to Fig. 1 (which shows the corresponding growth rates). 
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Fig. 3. — Dimensionless stress, T^^, of linearly unstable modes as a function of location, R, 
in the vicinity of the classical marginally-stable orbit (i?ms = for both magnetized and 

unmagnetized disks. The two groups of three curves represent stresses (i) for the fastest growing 
mode in a magnetized disk (top 3 curves) and (ii) for (wavenumber independent) unstable modes in 
a hydrodynamical disk (bottom 3 curves, interrupted at or around 3Rg). In each group, the solid line 
is for the case without radial stratification [N'^ = 0), the long-dashed line is for a stabilizing radial 
stratification (N'^ = O.IO^) and the short-dashed line is for a destabilizing radial stratification 
(A/"^ = —O.irP). The dotted line, which represents the stress for a long-wavelength mode with 
{kvA/^y = 10~^ in a magnetized disk, remains continuously positive throughout the region of 
interest even for this "hydro- like" mode. Note that, while signs are accurate, stress singularities 
for the bottom short- and long-dashed lines are only apparent: one can show that the dimensional 
stress vanishes, rather than diverges, at marginal stability (see text for details). 
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Fig. 4. (a) Dimensionless frequency, cj/Jl, of incrtio-gravity-Alfven waves as a function of 
dimensionless wavenumber, /cu^/O, in a magnetized disk without radial stratification. The 
various curves correspond to different locations in the vicinity of the classical marginally-stable 
orbit (same notation as Fig. 1). Long- wavelength inertio-gravity-like waves have nearly-zero (but 
strictly positive) frequencies inside of the classical marginally-stable orbit. The behavior of short- 
wavelength Alfven-like waves, on the other hand, is largely unaffected by location, (b) Dimensional 
frequency, w, of inertio-gravity-Alfven waves with specific wavenumbers as a function of location, R, 
in the vicinity of the classical marginally-stable orbit {Rms = 3i?g). From top to bottom, the various 
curves correspond to (fcfA)^/^^^ = 10~^, 3 x 10~^, 2 x 10~^, 10"^ and 10~^. The orbital frequency 
has been normalized to = 1 at i? = 3Rg in this plot and there is no radial stratification. Notice 
that, contrary to the purely hydrodynamical case, the position of the frequency maximum is not 
fixed at R = ARg but it varies with wavenumber; the frequency maximum disappears as magnetic 
tension becomes important. 



